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Abstract. In this paper we obtain some properties of orthogonal polynomials given by a weight 
function which is a limit density of a rescaled discrete-time quantum walk on the line, f f 

1 Introduction 

A quantum walk is the quantum analog of a classical random walk. Quantum walks are 
expected to play an important role in the field of quantum algorithms. A number of benefits 
for such walks are already known. Reviews and books on quantum walks are Kempe [6], 
Kendon [7], Konno [TOl [11], Venegas-Andraca [I5], for examples. There are two types of 
quantum walks. One is the discrete-time walk and the other is the continuous-time one. 
Here we focus on the case of a discrete-time walk on Z, where Z is the set of integers. 
Ambainis et al. |3] investigated the quantum walk intensively. In the present paper, we 
consider some properties of orthogonal polynomials given by a weight function which is a 
limit density of a rescaled discrete-time quantum walk on Z. Recently Cantero et al. [1] 
showed that the theory of matrix-valued orthogonal polynomials associated with a certain 
kind of unitary matrices, i.e., the CMV matrices, is a natural tool to study the discrete-time 
quantum walk on the line. It would be interesting to know a relation between their approach 
and our results. 

The rest of the paper is organized as follows. In Sect. 2, we define the quantum walk 
and explain the weak limit theorem. Section 3 treats a symmetric qunatum walk. In Sect. 
4, we give a result for an asymmetric case. Section 5 is devoted to a general case. 
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2 Quantum walk 



First we give a definition of one-dimensional discrete-time quantum walk. The time evolution 
of the quantum walk is defined by the following matrix: 



a b 
c d 



e U(2), 



where a,b,c,d G C and U(2) is the set of 2 x 2 unitary matrices. Here C is the set of 
complex numbers. The unitarity of U gives 



1, ab + cd — O, c — —Ab, d — Aa, 



where z is the complex conjugate oi z E C and A = ad — be. The quantum walk is a quantum 
version of the classical random walk with additional degree of freedom called chirahty. The 
chirality takes values left and right, and it means the direction of the motion of the particle. 
At each time step, if the particle has the left chirality, it moves one step to the left, and if 
it has the right chirality, it moves one step to the right. Let define 



1^) 



1 




\R) 
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so U acts on two chiralities as follows: 

U\L) ^a\L) + c\R), U\R) ^b\L) + d\R), 

where L and R refer to the right and left chirality state respectively. Here the set of initial 
qubit states is defined by 

^={^ = ^[a,P]eC':\a\' + \P\' = l}, 

where T is the transposed operator. Let Xn{= X!^) be the quantum walk at time n starting 
from the origin with the initial state e To explain more precisely, we introduce P 
and Q given by 



a b 




, Q 




c d 



with U — P + Q. The important point is that P (resp. Q) represents that the particle moves 
to the left (resp. right). Let m) denote the sum of all paths in the trajectory consisting 
of / steps left and m steps right. In fact, for time n — I + m and position x — —I + m, we 
have 

lj,mj 

summed over all non- negative integers Ij and rrij satisfying • — I and mi + - ■ •-\-mn — 
m with Ij -\-mj — 1. The definition gives 



S„+i(/, m) = P -l,m) + Q E;„(/, m - 1). 
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The probability that X„ = x is defined by 

PiXn = x) = ||S„(/,m)(^||2. 

We should remark that there is a strong structural similarity between quantum walks and 
correlated random walks, see Konno [12]. A typical example of the quantum walk is the 
Hadamard walk defined by the Hadamard gate U = H: 



V2 



1 1 
1 -1 



The walk has been extensively investigated in the study of the quantum walk. 

Quantum walks behave quite differently from classical random walks. For example, in 
the classical case, the probability distribution is a binomial one. On the other hand, the 
probability distribution of the quantum walk has a complicated and oscillatory form. For 
the classical case, the well-known central limit theorem holds. For the quantum case, a 
corresponding weak limit theorem was shown by Konno [HI E] in the following way: 

Theorem 2.1 If abed ^ 0, then 

lim P{u < Xn/n < v) = {1 — c(a, b, v?)^} k(x : \a\) dx, 

where 



kix ■ r) - _ /, Jx) da b-u^)- \a\' - \6\' + """^^ + """^^ 
i^yx . r ) — i(^^r,r)\X), cya, u . ip) — \a\ \p\ -t- 



7r(l — x^^\Jr'^ — ' |a 
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and if = Here Ia{x) = 1 (x G A), = (x ^ A). In other words, the limit density of 

Xn/n is foo{x) = {1 - c(a, b, (f )x} k{x : \a\). 



3 Symmetric case 

In this section we consider the symmetric case of the limit density foo{x), i.e., c(a, b:(f) = 0. 
Let /i be the probability measure on the real line M, with the density k{x : r), where < r < 1. 

Theorem 3.1 LetG^{z) := '^^^fj he the Stielties transform of n. Then 



, , zi^z"^ — r^) — -\/l — r'^\/ z^ — 

= (z2 _ 1)(^2 _^2) • 

Moreover, admits the following expansion as continued fraction: 



1 - VT^ 

z 



{^/T^-l+r^) /2 



rV4 



rV4 
z 
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Before the proof we will derive some consequences. 
Corollary 3.2 The monic orthogonal polynomials for ^ are given by: Po{x) = 1, 

xPnix) = P„+i(x) + 7„_iP„,_i(x), 

where 



2 



7o = 1 - Vl - r2, 7i = ^ , 1n = — forn> 2, 

under convention that P„i(x) = and 7„i = 0. 

In particular, for the Hadamard walk case (r = 1/ \/2) we have 

2- V2 V2-1 1 
70 = ^ , 71 = — ^ , ^" " 8 - 

Then we can compute a few first orthogonal polynomials: 

Po{x) = 1, Pi(x) = X, 

P2(x) =X' + , 

P3(a;) = + X, 

4 , -7 + 2V2 

P^{x} = x + X + , 

/ N 5 -4 + V2 3 7-3^2 
P^{x) =x^ + x' + — X, 



6 -9 + 2V2 4 21-8v/2 -2 + 
P6(x) = x H — X H — X + 



23 26 2^ 

From Corollary 13.21 we can easily compute the following generating function: 



Q{x,z)=Y,P^ 



00 

'"..[x]z'\ 

n=0 



Corollary 3.3 

(4 - Axz + r^z^) Q{x, z) =4 + (r^ -A + Wl-r^^ z"^ + (2 - 2Vl-r^ - r^) xz^. 
In particular, for the Hadamard walk 

(4 - 4xz + ^z^^ Q{x, z)=4+(^-^- + 2V2^ + [I ' ) 

The following fact is useful for computing G^{z). 
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Lemma 3.4 Define 



z 



rV4 
z 



Then we have 



2z 2 



A{z) = — 1/^2 _ 



^2 

Proof. The definition of ^4(2;) yields 

^^""^ " 2-rM(2)/4' 

so we liave 

r'^A{zY - 4:zA{z) +4 = 0. 

Tlierefore 

2z . 2 



wliere tlie sign is clioosen in sucli a way tliat A maps tlie upper lialf plane into the lower half 
plane, see [21 [S]- 

Proof of Theorern VJ . 1\ Denoting the continued fraction by G{z) we have 
G{z) = 



1 - 



{^/T^-l+r^) /2 
z - A{z)r'^/A 

1 



Vl — — 1 + 
z 

z + y z'^ — r'^ 

1 



z — 



z — {\/\ — r"^ — 1 + r^) [z — \/ z^ — r^) / 
+ Vl — r'^\J z"^ — 



^2 y^2 



2;2 _ 7-2 _j_ 2:1/1 — r^V 

z (z^ — r^) — Vl — r'^\fz^ 

(^2 _ 1) (^2 _ ^2) 
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which is equal to the right hand side of the first formula. Now, using the standard technique 
(see [21 E]) can verify that G{z) is Stielties transform for /x, which concludes the proof 
of Theorem 13.11 

Next we consider the mth moment of the measure /i: 



:= / x"'k{x:r)dx = / — dx. 

Jr j-r 7r(l - a;^) V?" - a; 

Theorem 3.5 Form >0 we have S2m+i(/^) = and 



The moment generating function is equal to 



1 — r'^z'^ — z'^y/l — r'^y/l — r'^z'^ 



(1- z^)(l-r^z^) 

We note that this result for the Hadamard walk case appeared in Konno et al. |13] . 
Proof. First we note that 



M,{z) = G,{l/z)/z = ^i_,2^^i_r2,2^ • 

Now denoting the sequence given in theorem by Sm, and its generating function by M{z), 
we have 

/2a fr'^"' 



^2n 



4 

2A;\ / 



n=Q L k=0 ^ ^ 

k=0 n=k+l ^ ' 



1 z^VT 



(l_^2) /l_4r^ 



1 — r^z^ — z'^y/1 — r^Vl — r'^z'^ 
= (l-z2)(i_^2^2) ^ 

which is equal to M^{z). We used the well known formula: 



E 

fc=0 



2a . 1 



x' 



k J a/1 - 4a; 
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4 Asymmetric case 

In this section, we consider an asymmetric case of the hmit density foo{x)- Denote by /i(r, c) 
the probabihty measure on M which has density (1 + cx) k{x : r), where r G (0,1) and 
c G [— 1/r, 1/r]. We note that c = leads to the symmetric case. 

Theorem 4.1 For the moment generating function and the Stielties transform of fi{r,c) we 
have 



(1 — r^2;^)(l + cz) — {z + c)z\J\ — r'^\J\ — r'^z^ 



{z"^ — r'^){z + c) — (1 + cz)^/l — r'^\J z^ — 



Proof. First we note that for the moments of /x(r, c) we have 



, , /■ vT^^r2(l + cx) , 

Sn[fi[r,cjj = I dx 

7r(l — x^) vr^ — x^ 



dx + c , dx 



-r 7r(l — x2)-\/r2 — x^ — x^)vr^^^^ 

s„(yu(r, 0)) + c ■ s„+i(/i(r, 0)). 



Hence we can use Theorem 13.51 to get 

oo 



n=0 



(1 — r^2^)(l + cz) — (-2 + c)zy/T^^r^y/T^^r^ 

= (1 -22)(1 _^2^2) • 

From this we obtain the Stieltjes transform G^(^r,c){z) = M^(^r,c)i)-I z) / z of ^{r,c). 
Denote by Pn{r, c), ■jn{r, c) the Jacobi coefficients of /i(r, c), so that 

1 

G';i(r,c)(^) = 



z - (jo{r, c) 



z - fji{r, c) 



^ . . 72 (r, c) 
z- P2{r, c) 



Since fi{r,—c) is the reflection of fi{r,c) we have Pn{r, —c) = — /?„,(r, c) and 7„(r, — c) = 
7„(r, c), so we can assume that c > 0. Using combinatorial relations between moments and 
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Jacobi coefficients (see [Il[l6]) one can check that 

/3o(r,c) = c(l - s), 



7o(r,c) = {l-s){l-c' + c^s), 

-c(l -s)(2-2c2-s + 2c2s) 



/?i(r,c 
7i (r,c 



2(1 - c2 + c^s) 
g(l-g)(2-2c^ + c^g + c^g^) 

4(1 - C2 + C2s)2 



where s := y/T—r^, and these coefficients are getting more and more comphcated. It is 
possible however to find them in particular cases. 

Proposition 4.2 



1 — -\/l — r^, n = 0, 

Pn{r,l)= { -{l-VT^)/2, n = l, 
0, n>2, 



ln{r, 1) 

Proof. We have 



vr^72(i - vr^72), n = o, 

rV4, n>l. 



z- {1- Vl - r2) 



v^r^72 (1 - x/F^) 



+ (1 - VT^^ /2 



r2/4 



r2/4 

z 



r2/4 
2; 



1 



;z + (1 - v^r^r^) /2 - A{z)r^/A 
1 



+ (1 - Vl^) + 



^2 ly 2 



+ (1 - v^r^) + V 



^2 .y^ 2 



2-2 _ 



^2 - r2 + (2 - (1 - v^r^^)) V- 
2 (1 - yr^) ((^2 - r2) - v/r^V-z' - r2) _ 
2 (1 - yr^) (z2 _ r2)(z ^T) ~ 



In a similar way one can prove 
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Proposition 4.3 



(1 - v^r^r2)/r, n = 0, 

/?„(r, 1/r) = <! -{l-^l^f/{2r), n = 1, 
0, n>2, 



7n(r, 1/r) 



vr^r2(i - V^^^)^ n = 0, 

rV4, n > 1. 



5 General case 



We consider a general case with the Jacobi coefficients {7n}^o "which are given by 

PO = To, Pi = 7l, ■ ■ ■ , Pn-l = 7n-l, P = 7n = 7n+l = • • • ■ 

The corresponding Stieltjes transform is denoted by G^'^\z). In this paper, we call the case 
(po,Pi, • • • ,Pn-i;P)-case. As wc showed in Sect. 2, the symmetric case induced by the 
quantum walk is a special case for n = 2, i.e., (1 — \/l — r"^, Vl — t2(1 — ^1 — r2)/2, r^/2)- 
case. In a similar fashion, we can obtain an explicit form of G*^"-* (2;) as follows. 

Theorem 5.1 



z 



where Uk{z) = zllk-i{z) - Pn-{k~i)^k-2{z) with Ilo{z) = z - Pn-iA{z),Il_i{z) = 1. i^ere 
/n particular, 

1 {2p-po)z -po\/z^ - Ap 



G('\z) 
G^'\z) 
G^'\z) 



{p-po)z^+pl 



{2p - pi)z +pi-\/z'^ - Ap 
{2p-pi)z'^ - 2pop + piZy/z^ - Ap 

{2p - p2)z- - -Ipip + p2r- \/rJ - ip 



{2p - p2)z^ + (P0P2 - 2poP - 2pip)2; + p2{z'^ - Po) V^j^ - Ap 

Applying the Stieltjes inverse formula, we have the following absolutely continuous part of 
the corresponding probability measure. 

Corollary 5.2 



(1)/ ^ 1 PQ^Ap-x'^ 

' ^'^^2^ ip-p,)x^+pl 



(2)/ N _ _± PoPl^yAp-x'^ 

^ ^''> 27r (p-pi)x4 + {po(pi-2p)+p?K + p2p^(-2v^'2^)^^^ 



(3)/ X _ i PoPiP2a/4p - x'^ 

P ~ 27r {p-p^)x^ + c,x^ + c2X^+pIpI ^(-2^.2^)^, 
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where 

Cl = (PO + Pl)(P2 - 2p) + (po + P2)P2, C2 = (po + Plfp - P0P2iP0 + Pi + 2^2)- 

For n = 1 case, see Obata [H]. In particular, when po = 2N and p = 2N — 1, p*^^^(a;) is a 
constant multiple of the density function of a Kesten distribution. 
Finally we consider the following asymmetric case: 

Po = 7o, Pi = 7i, P = 72 = 73 = ■ ■ ■ , qo = Po, q = Pi = p2 = ■■■ ■ 
In a similar way, we obtain 
Proposition 5.3 

Q{2,asym)^^^ ^ (2p - Pl)z - q{2p - pi) + Pi ^(^ - g)^ - 4p 

(2p-pi)(2; - qo){z - q) - 2pop + pi{z - go)\/(^ - q)^ - 4p' 

^(2,asym). ^ _ i PoPi^Ap- {x-qy 

27r (p -pi)(x - qo^ix - qY + {po(pi - 2p)(a; - q) +pf(x - qQ)]{x - go) +PoP 
X -^(9-205,(7+2^) (a;). 

We note that if go = = 0, then -;~) (^^esp. becomes G^'^\z) (resp. 
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